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ABSTRACT 

A  refinement  of  the  results  obtained  for  optimum  energy 
constrained  jamming  of  digital  receivers  is  obtained  by 
modeling  the  jammer  as  a  random  process.  In  the  modeling 
process,  a  random  time  arrival  or  random  frequency  errors 
are  accounted  for  by  including  these  effects  in  the  repre- 
sentation of  the  jamming  waveforms.  Performance  analyses 
are  carried  out  in  order  to  determine  the  effect  of  random 
time  of  arrival  and  random  frequency  errors  on  the  part  of 
the  jammer,  on  the  receiver  probability  of  error. 

The  mathematical  results  derived  are  programmed,  evalu- 
ated on  the  computer,  and  compared  against  ideal  optimum 
energy  constrained  jamming  strategies  previously  studied. 
Results  for  both  coherent  and  incoherent  receivers  are 
derived  and  analyzed  utilizing  conventional  binary  modula- 
tion schemes.  Results  show  that  generally  some  but  not  a 
great  deal  of  jammer  effectiveness  is  lost  due  to  random 
time  of  arrival  or  random  frequency  errors  associated  with 
the  jammer  waveform. 
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I.  INTRODUCTION 

The  structure  and  performance  of  digital  receivers 
operating  in  additive  White  Gaussian  Noise  (WGN)  is  well 
known.  While  this  noise  assumption  is  often  valid,  in  many 
cases,  especially  when  there  is  intentional  jamming  in  the 
channel,  the  WGN  interference  assumption  breaks  down.  For 
this  reason,  there  is  interest  in  determining  the 
performance  of  digital  receivers  designed  to  operate  in  a 
WGN  environment,  that  must  however  operate  in  a  different 
environment . 

In  this  thesis,  the  performance  of  digital  receivers  has 
been  analyzed  under  the  assumption  that  the  interference 
consists  of  additive  WGN  and  some  intentional  jammer 
waveform,  whose  model  represents  a  refinement  over  previous 
work  in  this  area.  It  is  to  be  expected  that  the  presence  of 
any  jammer  without  such  prior  knowledge  on  the  part  of  the 
receiver,  will  cause  the  performance  of  the  receiver  to  be 
degraded. 

For  the  case  in  which  the  jammer  waveform  is 
synchronized  with  the  digital  signal  and  has  exact  knowledge 
of  the  frequency  of  operation,  results  on  receiver 
performance  have  already  been  obtained  and  analyzed  [Ref 
1,2].  This  thesis  endeavors  to  investigate  the  effect  of  a 
more  realistic  jammer  model  in  which  the  jammer  lacks 
synchronism  with  the  digital  signals  as  well  as  exact  signal 
frequencies  knowledge.  This  model  is  used  in  conjunction 
with  both  coherent  and  incoherent  binary  digital  receivers 
designed  to  operate  in  a  WGN  environment. 

In  Chapter  2,  the  performance  of  coherent  receivers  is 
investigated  under  the  assumption  that  the  jammer  waveform 
lacks  synchronism  with  the  digital  signals  or  that  it  lacks 
exact   knowledge   of   the   frequencies   of   operation.   A 
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mathematical  model  of  the  jammer  waveform  is  introduced 
based  on  previous  results  on  optimum  energy  constrained 
jammer  waveforms  used  to  degrade  the  performance  of 
different  binary  modulation  receivers.  The  analysis  of  the 
performance  of  the  receivers  is  carried  out  and  the  effect 
of  miss-synchronization  (i.e,  timing  errors)  and  of 
frequency  offsets  (i.e,  frequency  errors)  are  discussed  in 
Chapter  5. 

In  Chapter  3,  the  work  carried  out  in  Chapter  2  is 
repeated  for  incoherent  receivers.  Specifically,  incoherent 
Amplitude  Shift  Keying  (ASK)  and  incoherent  Frequency  Shift 
Keying  (FSK)  receivers  are  analyzed,  performances  are 
evaluated  and  the  results  are  discussed  in  Chapter  5. 

In  Chapter  4,  graphical  results  are  presented 
corresponding  to  the  numerical  analyses  that  have  been 
performed.  In  most  cases,  the  graphs  display  receiver 
probability  of  error  as  a  function  of  signal  to  noise  ratio 
for  given  values  of  jammer  to  signal  ratio.  The  effect  of 
jammer  miss-synchronization  or  of  frequency  offsets  is  then 
analyzed  with  the  aid  of  the  graphs. 

Chapter  5  presents  conclusions  pertaining  to  what  was 
learned  from  the  results  obtained  and  the  numerical 
evaluations  carried  out. 
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II.  COHERENT  RECEIVERS 

A.   GENERAL 

The  correlator  receiver  structure  depicted  in  Fig  (5.1) 
which  can  be  shown  to  be  equivalent  to  the  optimum  single 
correlator  receiver  of  Fig  (5.2)  is  known  to  be  optimum  (in 
minimum  probability  of  error  sense)  for  discriminating 
between  two  signals  s-^(t)  and  SQ(t),  received  in  the 
presence  of  additive  White  Gaussian  Noise  (WGN),  with  0<t<T. 
We  define 

r(t)  =  received  signal  in  the  interval  0<t<T 

sd(t)  =  s1(t)  -s0(t) 

NQ/2  =  Power  spectral  density  level  of  the  WGN. 

P  =  probability  that  Sg(t)  is  transmitted. 

,  .  _2_  .   ln<-JL_>  ♦  i  (Ei-e0> 

2 

where  y  is  the  receiver  threshold  level  used  to  decide  on 
whether  Si(t)  or  SQ(t)  was  transmitted  and  E^  is  the  energy 
of  the  transmitted  signal  s^(t),  that  is 
rT 


e      =  r  s  <t>2  dt        i  =  0fl 

J0 


The  probability  of  error  P  of  this  receiver  is  derived  in 
many  text  books  [See  for  example  Ref  3] .  The  coherent 
digital  communication  receiver  of  Fig  (5.2)  can  be  analyzed 
in  terms  of  the  resulting  Pe  when  r(t)  contains  a  jammer 
waveform  n^(t)  in  addition  to  the  noise  and  Sg(t)  or  s-i(t). 
Thus,  under  these  conditions,  the  received  signal  appearing 
at  the  front  end  of  the  receiver  is  mathematically  described 
by 

r(t)  =  si(t)  +  n(t)  +  nj(t)        0<t<T   i  =  0,1 

12 


where  again  SgCt)  and  s^(t)  are  the  two  signals  used  to 
transmit  the  binary  information,  n(t)  is  a  sample  function 
of  a  White  Gaussian  Noise  process  having  a  power  spectral 
density  level  of  NQ/2  (Watts/Hz),  and  n^(t)  is  the  jammer 
waveform  present  during  the  signaling  interval  [0,T]. 

The  receiver  of  Figure  (5.2)  has  been  analyzed  in  so  far 
as  the  effect  of  n:(t)  on  the  receiver  Probability  of  error 
(Pe)  is  concerned  under  the  assumption  that  n^(t)  can  be 
modeled  as  deterministic  waveform.  [Ref  4] 

It  has  been  demonstrated  that  under  such  an  assumption, 
the  receiver  P   becomes 


P   =  P 

e 


erf  c 


*'   *  4llSJ|2  "   (ni'SH> 


I 


No 


♦  <1-P>  erf 


2    llSdH 


Ho 


I  I*- I  I 


(2.1) 


where 


<n.,Sd> 


"V*  ■  I] 


f 


n.(t)  •  S.(t)  dt 
J        d 


s.(t)   dt 
a 


With  P=l/2  ,  we  obtain 


■MI 

1  JS  fN 


«[V<] 


r  °  I* 

J    -co 


T~~2n" 


4    5n 


x 
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X 
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dx 


dx 


(2.2) 
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where 

f              2               1                          llSdH2 
S_    =      2  N      =   5 d   =    <n.    c    > 

q  I    No||Sd||2    J  d  2  J'    d 

B.   JAMMER  OPTIMIZATION 

In  the  previous  section,  Equation  (2.1)  allows  for  the 
evaluation  of  performance  of  the  receiver  under  the 
assumption  that  a  jammer  waveform  is  present.  From  the 
jammer  standpoint,  the  optimum  jammer  waveform  must  be 
chosen  in  such  a  way  that  it  maximizes  the  receiver 
Probability  of  error.  By  evaluating  the  first  derivative  and 
second  derivative  of  P  with  respect  to  'd'  which  is  the 
cross  correlation  between  the  jammer  waveform  and  the  signal 
difference  Sj(t),  the  optimum  jammer  waveform  can  be  chosen. 
Carrying  out  the  mathematical  operations,  we  have 

2,„    2.  .2 


S 

q 


[v*dl 


dP  S 


ad  T2n 


sinh(S         N.d) 
q         d 


3 


-  3„      _  Vlvfl 


d    P  S       N, 

e  q       a 


d    d2  T2iT 


>    © 


From  the  above  results,  it  can  be  seen  that  Pe  is  a 
monotonic  function  of  'd'  and  making  'd'  as  large  as 
possible  in  magnitude  results  in  the  largest  possible 
increase  in  P  .  In  the  limit,  as  |d|  —  «  we  have  Pe  —  1/2. 
However,  from  the  Cauchy- Schwarz  inequality 

I  |d|  |  =  |  |(nj,sd)|  |  <  |  |iij|  |  |  |sd|  |  (2.3) 

with  equality  if  n.(t)  is  proportional  to  s^(t).  Defining 

1/2 
I  I  n  j  |  |   =   PRj  '^   ,   where  Pnj   is   the  jammer  energy,   the 
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condition  |  d  |  -«.  «  implies  that  PRj  —  « 


when 


Under  the  condition  that  the  jammer  power  be  constrained  to 
P_  .•  ,  in  order  to  get  the  maximum  P  ,  Equation  (2.3)  can  be 
made  into  an  equality  by  setting 

rijCt)  =  K  sd(t) 

where  K  is  a  constant  of  proportionality.  Since  |  In^l  |  = 
Pnj  ,  K  must  be  set  to  the  value  PRj 1^2/ | | sd | | .  Thus  ||d|| 
is  maximized  by  setting 


n.(t)  = 
J 


1 11.1 1 


SJ<t) 
a 


(2.4) 


and   this   results   in 


being   maximized.   Substituting 


Equation  (2.4)  into  the  probability  of  error  expression  of 
Equation  (2.2)  for  P=l/2  yields 
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(2.5) 


By    defining 

E/NQ   =  SNR   ;  signal  to  noise  ratio. 

Pnj/E  =  JSR   ;  jammer  to  signal  ratio 
and  observing  that 

2   rT  9 

||Sd||       =  LS^t)    -    S0(t)3i£dt    =    2E(1-P) 
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where  p  is  the  cross  correlation  coefficient  between  two 
information  signals,  the  probability  of  error  expression  of 
Equation  (2.5)  becomes 


--  -  *  i  J 


2 
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V^NRC/T^p  -y2JSR) 
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-v^NRCyi^p  *y2JSR)        _x 

1      2 


dx 
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e      dx        (2.6) 
,/2F  1 


C.   EFFECT  OF  DETERMINISTIC  JAMMING  WAVEFORM 

In  the  previous  section,  the  design  of  the  jammer 
waveform  and  its  degrading  effects  on  the  receiver 
performance  has  been  analyzed.  The  effect  of  such  jamming 
waveforms  on  the  receiver  performance  has  been  found  in 
terms  of  receiver  probability  of  error.  From  Equation  (2.6), 
it  can  be  easily  seen  that  the  three  independant  variables 
which  affect  the  performance  of  the  receiver  are  SNR,JSR  and 
cross  correlation  coefficient  between  the  two  signals  used 
to  transmit  binary  information,  denote  by  p  In  a  binary 
signaling  set,  it  can  be  shown  that  -l<p<l.  For  Phase  Shift 
Keying  (PSK)  ,  p  =  -1  and  for  FSK,  p  =  0.  Thus  P£  for  PSK, 
using  Equation  (2.6)  is  given  by 


Hi  " 

J^§NR( 
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*   dx 
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V2-y2JSR) 
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-,/§NR(>/2*V2JSR>       — £— 

1      2       1 
e      dx  (2.7) 

-   CO 


and  for  FSK,  using  Equation  (2.6),  P   is  given  by 
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<n 
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D.   ANALYSIS  OF  THE  EFFECT  OF  THE  JAMMING  WAVEFORM  HAVING 
RANDOM  TIME  OF  ARRIVAL 

The  results  of  the  previous  section  demonstrate  the 
effect  of  an  energy  constrained  jammer  on  a  binary 
communication  receiver  using  a  strictly  deterministic  jammer 
model.  Such  a  model  however  must  be  refined  because  among 
other  imperfections,  the  jammer  signal  may  not  be  in 
synchronism  with  the  digital  transmission.  In  other  words, 
there  may  be  a  difference  in  the  time  of  arrival  of  the 
information  signal  and  the  jamming  signal.  The  modeling  of 
such  an  effect  will  be  accomplished  by  letting 

nj(t)  =  K  sd(t-T)    0^T-<T  (2.9) 

Observe  that  in  this  model,  the  optimum  jamming  waveshape 
has  been  maintained,  however  a  time  delay  parameter  f  has 
been  introduced.  In  order  to  be  able  to  analyze  how  the 
time  delay  7  influences  the  effectiveness  of  the  jammer,  and 
compare  the  resulting  receiver  Pe  with  that  involving  the 
idealized  jammer  model  in  which  the  jammer  operates 
synchronously  with  the  receiver,  we  will  require  that  K  be 
such  that 

llnil|2   ■  J   tK.S.(t-T)]2  dt    =P  (2.10) 

Therefore,  for  this  constraint  to  be  satisfied,  we  must  have 
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K  = 


nj 


Sj  (t-T)dt 
a 


(2.11) 


and  the  jammer  waveform  is  now  modeled  as 


n.(t)  =  K  S.(t-T) 
J         d 


(2.12) 


with  K  given  by  Equation  (2.11).  In  order  to  determine 
receiver  performance  using  this  new  jammer  model  we  may  use 
Equation  (2.1)  in  which  now 


(n.,S.)  b  K(S.    ,  3^) 
J   d        d,r  '   d 


(2.13) 


where 


<sdlT  •  V 


f 


Sd(t-T)  Sd(t)  dt 


while  the  remaining  terms  in  Equation  (2.1)  are  unchanged. 
Thus  the  Pe  for  the  receiver  operating  in  the  presence  of 
noise  and  a  (non-synchronized)  jammer  mathematically  modeled 
by  Equation  (2.12),  is  given  by 


P      =    P-erfc 

e 


+(1-P) «erf 


*•♦  4-iisji2-  k<s^  ,  .s,) 


d,T 


I 


No 


2         HSdll 

r    «'•  -g-  "5ill2-8(5d,T  >V    1 

I  I      Mr.  ....  J 


I 


Ho 


lis,  N 


(2.14) 


This  expression  can  be  evaluated  for  specific  modulation 
schemes,  such  as  Binary  Phase  Shift  Keying  (BPSK)  and  Binary 
Frequency  Shift  Keying  (BFSK). 
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1.   BPSK 

For  BPSK,  the  difference  signal  s^(t)  defined  on 
Page  12  becomes 

sd(t)  =  2A  cos<oct 

where  we  assume  that  0)CT  =  mr  and  fn'  is  an  integer.   Thus 


Sd,r 
so  that 


=  sH(t-r)  =  2A  cosio^Ct-r) 


f 


S  .  (t-T)dt 
a 


-f 


C  2Acos6>  (t-T)  l2    dt  *  2A2T  ■  4E 
c 


where 

A2T/2  =  E 

The  above  assumption  o;cT  =  n^-  has  been  used  in  order  to 
simplify  the  expression  "for  the  energy  of  Sj(  t-  t  ).  Thus 
from  Equation  (2.11)  and  Equation  (2.12), 


n.(t)  = 
J 


nj 


4E 


S.<t-T>  = 
a 


y/JSR    •SJ(t-T) 
a 


(2. 15) 


where  we  define 

Pn.:/E   =  JSR  =  jammer  to  signal  ratio 

In  order  to  numerically  evaluate  Equation  (2.14),  we  need  to 
specify  the  inner  product  term,  namely 


d.T    d     J, 


0 


2Acosu  (t-T)-2flcos«  t  dt 
c  c 


2A  T«cos6>  t 
c 


(2. 16) 
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By  substitution  into  Equation  (2.14),  receiver  Pe  for  some 
given  value  of  r  becomes 


P      =   P-erfc 

e 


+<1-P) -erf 


l 


lp  T^P  *  -B5  HsdH2-  4s  ''^'"V 

r^1  nsdii 

I|s„l1 


Ho 


<2. 17) 


where    furthermore,    since  p   =    -1» 
JL||Sd||2    =    2    SHR    (l-p>    =    4    SHR 


Jt£  llsdll2  -  2^m 

K    2A2T    -    4    SNR^TSR 


Ho 


Thus,  for  the  special  case  in  which  P=l/2, 


4SNR  -  4SNR^TSR'Cos«  t 

P  (t)  =  ^|erfc<  — — —  ) 

2,/2SNR 


■M' 


-4SNR    -4SNR,/JSRcosm    T 
2,/2SNR 


♦  erf<  —  )|       (2.18) 


Equation  (2.18)  can  be  expressed  in  terms  of  a  normalized 
delay  r^ •  Observe  that 

COS»  T   «  COSlu  T'-£—> 

c  c    T 

where 
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Thus  Equation  (2.18)  can  be  reexpressed  in  terms  of  r^ •  The 
integer  'nf  specifies  the  number  of  sine  wave  cycles  within 
the  observation  interval.  In  most  practical  cases,  n>10  , 
however  results  are  not  greatly  affected  by  this  parameter. 

It  is  apparent  that  in  most  cases,  tjj  cannot  be 
determined  a-priori,  so  it  must  be  treated  as  a  random 
quantity  with  some  associated  probability  density  function 
(p.d.f.).  A  logical,  yet  simple  choice  for  such  a  p.d.f.  is 
to  assume  t^j  uniformly  distributed  over  [0,1],  so  that 
Equation  (2.18)  must  be  interpreted  as  receiver  P 
conditioned  on  r^.  The  actual  (unconditional)  receiver  Pe  is 
obtained  by  integrating  over  the  p.d.f.  of  t^ >  namely 


P 

e 


■+£( 


4SMR    -    4SNR,/31jR.cos6>   tm 
erfc( £-*-> 


2y2SNR 

""  "->  1  "*. 


-4SNR    -    4SNR»/JSR.cos6>   t 

♦    erf( ;= ^L 

2>/2SNR 


BFSK 


For   BFSK,    we   have 
sd(t)    =    A    (cosio^t    -    COS(jQt) 
where   we    assume    that 
cosT    =    (w]_    +  0)q)T    =    2n7r 

and    'n'/m'    are    integers   with  n>m.      For   this    case, 

sd,T    =    sd(t_    t)    =    A    (coso)1(t-T)    -    coso>0(t-T)) 
so    that 

Sd    <t_T>    *t    -   A    I       Ccos»1(t-T)    -    cos«0(t-T)]2   dt 
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=  fl2T   «  2E  (2.20) 

From  Equations  (2.11),  (2.12)  and  (2.20)  the  jammer  waveform 
becomes 


J       4  2E 


Sd<t~r>  (2.21) 


Furthermore,    the    inner  product    term   takes    on   the    form 


J 


,S.)      =  fl2£COS«f  (t-T)      -     C05«ft(t-TH 

d,T'    d  J_  1  0 


[cosu.t  -  costt.t]  dt 

1  SJ 

2 
AT 

— = [costf   t    -    cos«   t]  (2.22) 


so     that     substitution    of    these     results     in    Equation     (2.14), 
yields 


In-: — =r+2SNR    -/2JSR-SNR*  (cosu.t-cosw^t) 


P    (t)    =    P-erfc 

2ySNR 


1-P 1 0_ 

2,/§NR  ' 

ln.P^-2SNR    -y2JSR-SHR. (cos«,t-cosm„t) 

2*/SHR 


L_L . 1 


(2.23) 

Utilizing       the       previously       defined       normalized       delay    r-^ , 
Equation    (2.23)    becomes 


ln^— =r+2SNR-y2JSR  •  SNR  •  (  2cosnrrr„  •  cosmnTw  ) 


P    (t    )=P-erfc 

2,/l3NR 


1-P H N_     | 

2^NR  ' 

p 

Iln-5 — 5-2SNR-y2JSR'SNR(2cosnnTu-cosiwiT    )     •, 
2^MR  * 

(2.24) 

where 
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0  <  Tfl  <  1 

Assuming  once  again  that  tm  must  be  treated  as  a 
random  quantity  with  uniform  p.d.f.  over  [0,1],  for  the 
special  case  of  P=l/2,  we  obtain 


VV-J,  W  dTH 


J0 


2SNR  -  2SNRy2 JSR  •  cosnnr  •  cosmrrr 

rfc.( 5 ?i 

2v^NR 


-2SNR    -    2SNRi/2JSR  •  cosnnr    •  cosmnT  , 

+erf.( ** *  )  dT„ 

2SNR  '   N 

(2.25) 

E.   ANALYSIS  OF  THE  EFFECT  OF  THE  JAMMING  WAVEFORM  HAVING 
RANDOM  FREQUENCY  ERRORS 

In  actual  situations,  it  is  oftentimes  difficult  for  the 
communication  jammer  to  know  the  exact  transmission 
frequency  or  frequencies  of  the  adversary's  communication 
system.  This  lack  of  knowledge  may  be  due  to  the  use  of 
imprecise  frequency  estimation  methods,  or  due  to  frequent 
frequency  changes  made  by  the  communication  system's  users. 
Thus  in  this  section  the  effect  of  lack  of  precise 
transmission  frequency  knowledge  on  the  part  of  jammer  is 
studied  and  analyzed  insofar  as  coherent  binary  digital 
receivers  is  concerned.  Specifically,  BPSK  and  BFSK 
modulation  only,  will  be  considered. 

1.   BPSK 

The  lack  of  precise  transmission  frequency  knowledge 
on  the  part  of  the  jammer  is  modeled  by  introducing  a 
frequency  difference  Ao>  in  the  known  optimum  jammer  model 
for  BPSK.  That  is , 

nj(t)  =  K  cos(Wc  +  Aa>)t  (2.26) 
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where  Aco  is,  in  a  sense,  the  jammer  frequency  error.   Since 
the  jammer  power  is  constrained  to  be  Pn-j»  we  must  have 


Hn.B* 

j  Au 


K2Tf 


[Bcos(»  +Au)tl  dt 

c 


sin2<«  *A*)T 


1  + 


\w  -rum  I   1    1 

__E 

<w  +  /W)T 


(2.27) 


Thus  we  have  from  Equation  (2.11), 


E  = 


2  P 


nj 


I '[ 


1  ♦ 


sin2<w  +Au)7 
c 


2(6*  +At* 

c 


>T  J 


(2.28) 


In  order  to  determine  the  effect  of  this  jammer  on  the 
coherent  receiver,  all  that  is  needed  is  the  determination 
of  the  inner  product  (n; , s^) a  .  Thus 


f 


(n..S.).       «     I       Hcos<«    +Au) t • 2Acosw    t«dt 
j'    d    At*         J0  c  c 

sin(2»    +Ao)T 
sKAT|sinA«T    + 


AuT 


(2o»  +  Au)T 
c 


(2.29) 


Defining  now  a  normalized  frequency  error,  Acom   where 


N, 


the  constant  K  and  (n-j>sd)AM  can  be  rewritten  as 


E 


2  P 


nj 


T  1  ♦ 


sm2w  T(  H  u„  ) 
c Nl 

2«  T(l**»„) 
c      N 


2 


and 
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'VVa.  "  KAT 


fl/T 


[ 


sin*  T«A«M    sin«  T(2*A*»  ) 
c     "4.      c       w 


o     N 
2  P 


nj 


«  T(2+46»„) 
c       N 

1_ 

2 


sin2<*  T(l+6»„) 


T  1  ♦ 


sin*>  T-A«M    sinw  T(2+A«u> 
c     N        c       N 


c     N 


6*  T(2+4«M> 

C         N 


(2.30) 


Substituting  Equation  (2.30)  in  Equation  (2.1),  we  obtain 

ln_P_  ♦  4SNR  -4SKHy3SRCB3CC3  , 
P  (A«)  =  P-erfc   J-^ZZZT" J 


♦  (1-P)«erf 


2y2SMR 
ln-y^  -  4SNR  -4SNRV^JSR-CB3-CC3 


2,/2SNR 


(2.31) 


where 


B  = 


2 


1  + 


sin«  Td+Au.,) 
c       N 


«  T(l+4u„) 
c       N 

sini*  T«4«„     sine*  T(2+4«„) 
c     N         c       N 


c     N 


6>  T(2+A«„) 
c       N 


Equation  (2.31)  yields  the  receiver  Pe  when  Ao>^  has  some 
known  value.  Unfortunately,  in  many  cases  Ao>^  will  not  be 
known  exactly  and  must  therefore  be  modeled  as  a  random 
variable.  Thus,  in  order  to  compute  the  average  value  of  Pe, 
it  is  assumed  that  the  probability  density  of  Aw^  is 
uniform  on  (Acot  ,Aco^),  where  A(0l  ^-s  t^ie  lower  limit 
and  A<x>h  is  the  upper  limit  of  the  normalized  frequency 
error.  With  this  assumption,  the  average  Pe  can  be 
calculated  as 
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e    **H-A»L 


A* 

H 

PeiAuHy'dA*H  (2.32) 

or  equivalently 

4SNR  -4SMR^JSRrB3  CC3 


P  = 

e  Aw 


— !-s C  P-erfo-   — 


2y2SHR 


In  P    -  4SNR  -4SNR*<JSRCB3 • CC1 
♦  (l-P)-erf 


]>d^N 


2y2SNR 

(2.33) 

where  B  and  C  have  been  previously  specified  and  are 
functions  of  Aco^j.  Results  on  the  evaluation  of  Equation 
(2.33)  will  be  presented  in  Chapter  4,  under  the  assumption 
that  P=l/2. 

2.   BFSK 

The  lack  of  precise  knowledge  of  the  transmission 
frequencies  on  the  part  of  the  jammer,  is  modeled  by 
introducing  frequency  differences,  Ac^i  and  Acoq  in  the 
known  optimum  jammer  model  for  BFSK.  That  is 

Ilj(t)     =     K[C0S(GJ;L     +     AcO]_)t     -     COS(COQ     +     AwqH] 

where  in  a  sense,  Acoi  and  Acoq  are  the  jammer  frequency 
errors.  Since  the  jammer  power  is  constrained  to  be  P  ^ ,  we 
have 

lln.ll2    =    P         =  K[cos(«.+4u,)t    -    cos(«/x+A«/a)t]2d1 

J  nj         J_l  11  0         © 


-_2_       r          sin2(«,+/iu1  )T  sin2  («„+4«_  )T 

2+    - - ♦    - - 

2<«i+/to,)T  2(6>.+^6>^)T 

11  V>          so 


sinU,  +Au,  ♦«_  ■*•£«_)  T              sin(«.  ■*•£«.  -«    -A6»_  )T 
1 10  0  _ 1 1 0  0 


( &» ,  +A<a .  ♦« _  +4«„ )  T  ( « .  ♦  A«* .  -«    -dw^ ) T 

1  1         id  V>  110V 


Thus 
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,.r — : ^ p 

^T(2+sinc2<«*,  +Au.  )T*sinc2(»^*4«_  )T-sinc(«    +4«*    >T-2sinc<«    -»  . )  )  J 
11  0         0  ss  da 


I2P 
nj 

•CD! 

J       T 

where 

sine    x    a 

sin 

X 

X 

cos   =   o>i 

+   <o0 

6>d    =    (0! 

+     O)0 

and 

D= 


T(2  +  sinc2<», +^«,  )T+sinc2(«_+A«_)T-sinc<»  +A<*    )T-2sinc  («*.-».))  J 
11  0    0  ss  dd 


In  order  to  determine  the  effect  of  this  jammer  on  the 
coherent  receiver,  all  that  is  needed  is  the  determination 
of  the  inner  product  (nj,s^).   Thus 

rT 

(n  .,S    )=      K£cos(»   +Att    >t-cos<4»0+A«»0>t3  •  A£cos«1t-cos«0t 3dt 
RflT         ,    sinA^T         sin(2«1+il«1)T         sin(«d-A«0)T 


2~~    '     I           A^JT  *            <2«»1»A*»1>T                 (6»d-A«0)T 

sin(«    +4«*->T  sin(«.*A«1  )T  sin(«    +£«,  >T 

5         0  d         l  s i 

(«    +Au^)T  (u.+iW.  )T                 («^+^«.)T 

s         0  d          l                             s          l 

sin<2«    +ite_>T  sinA«_T 

0         0  0 

*      TT^A^yf  AvJ 


We  define  now  two  normalized  frequency  estimation  errors, 
namely 

A(A 

*»»  •   -ST"     i=0'1 

1 

so  that 

Aco]_T   =   o>]_T  AwN    ]_ 
AwqT   =    coqT  Aco^,0 
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and 

<W    "   -T1   •    £"  (2.34) 

where 

F    =    5inc(«,T -4«u    t)    +    sinc(»fT- (2*^«M    .)) 
~    sinc<«dT-«0T   ^M>0)    -    sinc(«sT-«0T.iWMj0) 

at  H,  1  si  N, 1 

♦    sincCttgT-itejj   0>    +    sinc<«0T<2+4«         )) 

Using  again         Equation  (2.1),  Pe(Aa>N    i»Ao>jj    q)  for 

known   Aco^    i    and   Ao>n   q   values    become 

ln.J|=.  +2SNR    -2SNR^JSR-CD3 -CF3 

rlnj^=-  -2SNR    -2SNR^/TSR-CD3 -CF3, 

♦  (l-P)-erf    — ■ (2.35) 

1  2^HR  J 


If  we  now  treat  Aco^  i  and  Acojj  q  as  random  variables 
uniformly  distributed  over  appropriate  ranges,  and 
furthermore,  assume  that  the  normalized  frequency  errors  are 
statistically  independent,  the  probability  density  function 
of  Ao>n  i  and  Aco^j  q   becomes 

"'^h.i'^nV  =  <^m,i,h-^h,i,l>^',h,9,h-^k.o>l' 

where  AwN)x    is    in    the    region    (Aa)N}i?L  ,  Aa>N>]_jH), 

and  AcoN)q  is  in  the  region  (AcoN>o,L  '  A<°N,0,h)-  The 
average  P   is  finally  obtained  from 

*WNfl,H  *WNf0,H 


e  e    N,lp   N,0  r         Nfl*   M,0      N,l    N,0 
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or   equivalently 


"N.l.H^N,©,!! 


!=,iS-M,lIH-d-M(l,L)<4-H>0,H-^Mie,L>     L  J 


4**H,1,L   4"H,«,L 


P  «erfc 


ln-r^p    +2SNR    -2SHRv^TSR.  CD3 -IF3 


[__I 


2^§NR 


♦<1-P) -erf 


ln-r—-   -2SNR    -2SNRv^JSR- CD]  •  CF3 


[__I 


2VSHR 


l|dilttMlId^MfO 


(2.36) 


The  results  of  the  evaluation  of  this  Probability  of  error 
expression  will  be  shown  in  Chapter  4. 
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III.  INCOHERENT  RECEIVERS 

A.  GENERAL 

In  incoherent  systems,  the  phase  of  the  carrier  signals 
is  not  available  at  the  receiver  so  that  the  phase  must  be 
treated  as  a  random  variable  which  is  typically  assumed  to 
be  uniformly  distributed  over  [0,  2tt]  .  As  such,  we  may 
expect  the  performance  of  an  incoherent  receiver  to  be 
degraded  in  comparison  to  the  performance  of  the 
corresponding  coherent  receiver.  However,  because  of  their 
simplicity,  incoherent  systems  are  widely  used  in  many 
applications . 

The  analysis  to  be  carried  out  on  the  performance  of  an 
incoherent  receiver  in  the  presence  of  WGN  and  jamming  is 
based  on  previous  work  in  which  it  was  assumed  that  a 
deterministic  jammer  model  was  adequate,  and  that  the 
optimum  energy  constrained  jamming  waveform  for  the 
corresponding  coherent  receiver  can  act  as  a  good  jammer  for 
the  incoherent  receiver  also.  Thus,  such  near  optimum  jammer 
signals  are  studied  and  evaluated  in  terms  of  their  effect 
on  the  performance  of  incoherent  receivers,  under  more 
realistic  conditions  now  in  which  the  jammer  is 
miss- synchronized ,  or  lacks  exact  knowledge  of  the  signal's 
operating  frequencies. 

B.  EFFECTS  OF  OPTIMUM  JAMMING  WAVEFORMS 

1.   ASK 

Analysis  of  the  incoherent  receiver  starts  from  the 
assumption  that  r(t),  the  signal  appearing  at  the  front  end 
of  the  receiver  can  be  mathematically  modeled  by  either 

r(t)  =  A  sin(ajct  +  0)  +n(t)  +  nj(t)    0<t<T 

or 
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r(t)  =  n(t)  +  nj(t)  0<t<T 

where  0  is  a  random  variable  uniformly  distributed  over 
[0,27r],  n(t)  is  a  sample  function  of  a  White  Gaussian  Noise 
process  having  power  spectral  density  level  NQ/2  (Watts/Hz), 
and  n^(t)  is  the  jamming  waveform  present  during  the 
interval  [0,T] . 

In  the  absence  of  n^(t),  the  optimum  receiver  for 
the  binary  ASK  problem  being  presented  is  well-known.  Its 
derivation  is  well  documented  in  the  statistical  detection 
theory  literature  [Ref  3]  .  The  receiver  structure  is  shown 
in  Figure  (5.3).  An  equivalent  form  of  this  structure  is 
shown  in  Figure  (5.4). 

In  this  section,  the  effect  of  nj(t)  on  this 
receiver  is  analyzed  by  evaluating  the  resulting  P  ,  under 
the  assumption  that  n^(t)  is  a  deterministic  waveform, 
however  unknown  to  the  receiver  itself.  This  analysis  has 
been  carried  out  in  [Ref  2]  ,  but  is  repeated  here  for 
completeness  and  the  general  result  is  applied  to  the 
specific  problem  analyzed  in  this  thesis.  Receiver 
performance  evaluation  requires  determination  of  the 
statistics  of  either  G  or  G,  where  G^  is  the  output  of  the 
quadrature  detector  and  is  given  by 

G2  =  X2  +  Y2 

where 

X  ■  I   r(t)-sinw  t-dt  ■  (r,S) 
.  J0  c 

and 

Y  »     r(t)'00s*>  t-dt  m    (r,C) 

Provided  that  the  random  variable  0  is  fixed  to  some 
value  0  ,  X  and  Y  are  conditional  Gaussian  random  variables 
with 
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ECX|H1,«3  =  (flS^,S)  ♦  <nj»S)  =  ^xl 
EtY|Hlt^3  =  (AS^,S)  ♦  (nj»S)  =  "Vj 


a 

a  «■» 


where  Sn  is  used  in  place  of  sin(coct  +  Q)  .  Also,  it  can  be 
shown  that 

Var<X|Hlf0>  =  EC(n,S>2>  =  ^[i~ ^TT        =  VarCXlH0> 

and  similarly 

f    sin2«  T  , 

Var^Y|Hlf5>  =  =21  i-  —^    =  VarCY|H0> 

1         c   J 

For  convenience,  we  assume  cocT=n7r  ,  where  'n'  is  an 
integer.  Thus  the  'sine'  terms  above  vanish  resulting  in 
the  simplification 

VarCX|Hlf0>  =  Var{Y|Hlf0>  =  !!HZ  «  02 

If  this  assumption  is  not  made,  an  additional  term  remains. 
However  ,  if  (tt/coc)<<T,  the  additional  'sine'  term  is 
small  and  can  be  neglected.  Furthermore,  the  covariance 

Et  EX  -  ECX|H. ,0>3 «CY  -  E{Y|H.,0>]  |  H.,0> 

NoTf  A"  si"2"cT  I 
=  — 1     2.cT     j  "• 

provided  the  assumption  on  to  holds.  This  implies  that  for 
any  given  value  of  Q  ,  both  X  and  Y  are  uncorrelated 
Gaussian  random  variables  and  therefore  statistically 
independent.  The  density  function  of  G  conditioned  on  the 
phase  0  is  non-central  Chi-Squared  distributed,  and  is  given 
by 

g+M 
P     (g|H  0)    -  —j  ■•   2°   .1  {M   >  g>0 


G    \a  x       2a 


a 


2 
and  zero  otherwise,  where  <j   is  defined  above  and 
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M  =  E2CX|Hlf^3  ♦  E2CY|Hlf^3 

=  £  (AS^.S)  +  (n.,S)  J2  +  C  (ASflfC)  +  (n.,C>  J2 


Due  to  our  assumption  <drt  =  n7r  >  we  obtain 


AT  2 
M  *  <~>   ♦ATC<n.fS>cos0  ♦  <n.,C)sin0J 

♦<n.f S)2  ♦  <n.,C)2  (3.2) 

J  J 

If  r(t)  consists  of  noise  and  jammer  only,  then  X  and  Y  are 
also  independent  Gaussian  random  variables  with 

E<X|H03  =  (njfS)  =  i^ 
ECY|H03  =  (njfC)  =  ^ 

so  that  the  density  function  of  G   assuming  no  signal  is 
sent  is 

P    <g|V  .    *    .    *°2  .io<^> 

G  2a  a 


and  zero  otherwise,  where 

W»  *  (n  ,S)2  *<n   C>2  (3.3) 

J         j 

If  we  now  let 


P(H1>  =  Pr  (sinusoid  transmitted} 
P(Hg)  m    Pr  (no  signal  transmitted} 


then 
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Ht>  r 

J0 


Pe=  P<H1)    p  2(g|H1)dg  ♦  P(H0)     p  2<g|H0)dg      (3.4) 


vf 


where 


r  *° 

p  2<g|H.)  =     p   <g|H  ,0>.p  (0)  dtf 


i=0,  I 


oo  G 


The  second  integral  in  Equation  (3.4)  can  be  expressed  as 
follows 


D 


g+w 

2 


oo 


2a 


V2+M» 


o     2      * 
a 


V-e         I  (^i-  V)  dV  =  Q< 

o   a  a  *  a 


(  ^,-2-, 


2 

a 


where 


Q<a 


,0)  =  [V. 


2   2 
V  *a 


I  <aV)  dV 
a 


is  a  well-known  tabulated  function  called  the  Marcum  Q 
function.  The  first  integral  in  Equation  (3.4)  can  be 
expressed  as  follows 


2 

[      P  2<g|H1§0>  p^m-dfl  |-dg 

J0    **  —  oo 

2 
p  oo  -q 

P(a(o>)        p  2(g|HiV0>    d* 
J  -co  •*  0 


(3.5) 


where  the  Equation  (3.5)  has  been  obtained  by  changing  the 
order  of  integration.  The  inner  integral  of  Equation  (3.5) 
can  be  expressed  in  terms  of  the  Marcum  Q  function  as 
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f 

J0 


g+M 

2°2.I  (  £*->. 


2a' 


o  2 

a 


dg 


-J 

J0 


2 


V 


I  (_.V)  dV 
o  a 


l-Q(^5_|-iL_) 


Therefore   the   receiver  probability   of   error  Pe   can  be 
written  as 


'  =P(H  >  Q(^C-f-2— >  +P(Hf)| 
e     id      ^      a    '    a  1 


p2n 
J0 


^ 


>d# 


(3.6) 


where  the  dependence  on  0     is  imbedded  in  the  terms  M  and 

the  M'  ,  the  threshold  77  which  the  receiver  sets  assuming 

that  no  jammer  is  present  is  obtained  as  the  solution  to  the 
equation 


2 
AT 

No 


I  ( 
o 


2Ar? 
No 


>  =  X 


0 


0"   P(H.) 


which  can  be  equivalently  put  in  the  form 


A2T 
No 


I  ( 

o 


A2T 


)  =  X 


4    No 


0 


The  I(»  )  function  used  here  and  also  previously  used  in 
conjunction  with  the  development  of  Pe  is  the  modified 
Bessel  function  of  the  first  kind.  If  the  definition  of 
average  signal  energy  previously  introduced  is  used,  we  have 

E  =  (A2T)/4 
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which  is  reduced  in  half  in  comparison  to  the  signal 
transmission  case,  due  to  the  fact  that  the  information 
bearing  signal (s)  do  not  have  equal  energy. 

In  order  to  afford  comparisons  with  the  coherent 
receiver  case,  we  implicitly  boost  the  value  of  the  signal 
amplitude  A,  to  obtain 

E  =  (A2T)/2 

in  order  to  have  agreement  with  previous  cases  insofar  as 
signal  energy  is  concerned.  Thus  the  threshold  determination 
equation  now  becomes 

-SNR         

e      -I  (  y2SNR  --2-  )  =  X_ 
o  a  0 

If  we  assume  P(H1)  =  P(HQ)  =  1/2,  then  . 

and  the  threshold  setting  equation  becomes 

-SNR 

e      -I  (  ,/2  SNR  --2—  )  ■  1 
o  a 

If  we  are  to  find  the  optimum  jammer  waveform  so  as  to 
maximize  P  ,  an  attempt  must  be  made  to  solve 

dP                                                                 dP 
—  =  0         and  —  =  0 


Unfortunately  the  resultant  equations  are 
mathematically  involved  and  do  not  appear  readily  solvable 
for  n.:(t).  It  seems  however  that  a  good  jammer  waveform  can 
be  postulated  based  on  the  results  obtained  for  coherent 
ASK.  It  was  found  for  such  a  case  that  the  optimum  n^(t) 
under  the  constraint  that  the  energy  of  n^(t)  be  limited  to 
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some  value  P_^  is  a  tone  at  the  carrier  frequency.  Thus  the 
following  jammer  waveform  can  be  used  as  a  potential  near 
optimum  jammer,  namely 


n.(t)  = 
J 


nj  T 


sinu  t 
c 


0<t*T 


(3.8) 


Observe  that  with  this  choice,  ||nj(t)||  =  Pnj .  The 
probability  of  error  Pe  can  now  be  determined  using  the 
threshold  setting  equation  and  the  previously  derived 
expressions  for  M  and  M'  .  The  effect  of  the  near  optimum 
jammer  waveform  on  the  receiver  (i.e,  incoherent  receiver 
performance)  can  be  analyzed  by  evaluating  Pe  as  a  function 
of  JSR  using  Equation  (3.8). 

It  can  be  shown  that  when  P(Hq)  =  P(H1)  =  1/2  the 
probability  of  receiver  error  given  by  Equation  (3.6) 
becomes 

P   =  -4r~         H-Q<y2-SNR.JSR,-^-> 
e     2   [  '    a 

p2n   

-  -2—     Q(i    2-SNR(l  +  2,/JSR.cos0+JSR)  , -^-)  -d0    (3.9) 

a 


<b 


where 

JSR  =  Pnj/E 

In  section  C  and  D  of  this  chapter,  Pe  will  be  evaluated 
once  again  under  a  more  realistic  jammer  model  that  includes 
miss-synchronization  and  frequency  offsets. 

2.   BFSK 

For  binary  incoherent  FSK  with  a  jammer  present,  the 
received  signals  under  the  two  hypotheses  are  either 


H   ;  r(t)  =  A-sin(«.t  +  ®)  +  n(t)  +  n. 
1  .  *  J 


(t) 


0<t<T 


or 
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H^  ;  r(t)  =  A-sin(u0t  ♦  <$>)  ♦  n(t)  +  n.(t)        ©<t*T 


By  separating  the  frequencies  o)]_  and  <og  sufficiently,  we 
can  form  signals  that  are  orthogonal,  have  equal  energy,  and 
have  the  same  advantage  of  ease  of  generation.  Here,  Q 
and  ft  are  assumed  to  be  independent  random  variables 
uniformly  distributed  over  [0,2«jt]. 

The  receiver  function  is  to  compare  the  envelopes  at 
the  output  of  each  channel  once  every  T  seconds  and  decide 
in  favor  of  the  larger  of  the  two  envelopes  (Figure  5.5). 
It  is  assumed  that  the  probability  of  sending  either  one  of 
the  two  signals  is  0.5.  For  the  purpose  of  analysis,  assume 
first  that  a  'mark'  signal  has  been  transmitted,  that  is, 
the  hypothesis  Hi  is  true.  An  error  is  committed  if  qg 
exceeds  q-^.  An  error  is  also  committed  if  q^  is  larger  than 
qg  when  a  'space'  signal  has  been  transmitted,  that  is,  when 
the  hypothesis  Hg  is  true.[Ref  3] 

Let  PeT  denote  the  probability  of  the  first  type  of 
error  described  above,  which  is  expressed  as  Pr(qg>q-^  |  H-i). 
Under  the  assumption  that  a  'mark'  signal  has  been  sent,  the 
output  q-t    of  one  of  the  envelope  detectors  is  given  by 


*1 

where 


2  =  X-,2  +  Y-,2 


X   =     r(t)  sin«  t  dt         Y   =     r<t>  cosu  t  dt 

1    J0  l  J0 

Observe  that  X-i  and  Y-i  conditioned  on  the  phase  and  either 
of  the  two  hypotheses  are  Gaussian  random  variables  with 

EtX.lH.,03  =   Asin(uf t+0)  sinu.t  dt  +    n . (t ) • sino, t  dt 
S  CAS    fSl>  ♦  Cn.,8,) 
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and 


ECYf|Hlf0>  =   A-sin(«lt+5)'Costf,t'dt  ♦    n.(t)-cosw,t  dt 


5  <ASlf,fCl>  ♦  (n.,Cl) 


where  AS^  a  represents  the  function  Asin(co-j_t  +  0)  .  S^ 
represents  the  function  sinco^t  and  C^  represents  the 
function  costoit  .  Likewise,  assuming  again  that  o>^T  =  n-jr  , 
where  'n'  is  an  integer,  and  that  n(t)  is  a  zero  mean  white 
Gaussian  noise  process  with  P.S.D.  level  NQ/2  Watts/Hz,  we 
obtain 

VarCX1|H1,«3  «  Var*^  |HX ,0)  -   N°T 

Furthermore,  it  can  be  shown  that 

EC  [X   -ECX1 |Hlf^>3 .CYX  -ECYj JH^  , 03 3  |  H  ,8    3=0 

so  that  the  conditionally  Gaussian  random  .variables  Xi  and 
Yi  are  uncorrelated  and  therefore  independent.  The  sum 
involving  random  variables  X-j_  and  Y-j_,  and  producing  <3i  s 
will  result  in  a  non-central  Chi-Squared  distribution  so 
that 


P  2i      <q JH  ,0)    =  — J_.e    2*    .i  <_L_Li_)      q  >0 
Ql'Hl,*  2a2  °     a2  X" 


where 

.2,„  ...   ..    „2 


m 


1X  =  ExtX1|H1,^>  ♦  E^YjH^*! 

=  C(AS,  ofS.)  ♦  (n.,Sf)l2  ♦  C(AS,  olCt)  ♦  <n.,C,>32 
1 ,0*  i        J   I  1,0*1        J   i 
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and  o  -         N0T/4.    Using    standard    random    variable 

transformation  techniques,  it  can  be  shown  that 


ql+"il 


**i        or,2  q/  ™,  , 

QjH^ff  l       l  a2  a2  l 

so  that 


p2n 
PQ  <qjH  )  =     p        (q,|H,fl).p  C0)  dB 

1        J*  QJh,.* 


2 
0    a 


ql  *mll 
2n 


2  a2    T   'l^ll 


I  < 5 >-d0    (3.10) 

o      2 


where  the  dependence  on  0  is  imbedded  in  the  term  m-^.  On 
the  other  hand,  the  output  qg  of  the  other  envelope  detector 
(Figure  5.5)  when  Hn  is  assumed  to.  be  the  true  hypothesis, 
is  given  by 

q02  =  xo2  +  V 

where 

X0  *     r(t>  sin«0t  dt         Y^  =     r(t)  cos«0t  dt 

Following  a  similar  procedure  to  the  one  used  above,  it  can 

be  shown  that 

rT  fT 

E<X0|Hlf0>    =      ft-sint^t+d)  .sin«0t  -dt    ♦         n  .  ( t  )  •  sin«0t    dt 

S    (flSl,*'V    +    (nj'V 

and 
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E<Y/jH,.0>    o[  A-sin(«,t*«) -cos»_t  «dt    ♦         n.(t)«cos«_t    dt 
•'    l  J©  l  0  J©    J  0 

S    <AS1,*'V    +    <VV 

Here,  Sq  represents  the  function  sina>gt  ancl  Cq  represents 
the  function  cosont    .  It  can  be  demonstrated  that 

Var<Xg|Hlf0>  =  VarCY0 |HX ,«>  »  — ^ 

and  also  that 

so  that  the  conditionally  Gaussian  random  variables  Xq  and 
Yq  are  uncorrelated,  hence  independent.  Thus,  similar  to 
Equation  (3.10),  the  expression  for  the  conditional  density 
function  of  qQ  becomes 


*  ^iv-  4 


q«2*m»> 


0  2n 

v0 


q0      2  a2.  .  V  mOi 


^  — e         'I  ( 0    )-dg   (3.11) 

a  a 


where 

and  now  the  dependence  on  0    is  imbedded  in  the  term  mil, 

In  order  to  compute  the  probability  of  error,  we  can 
use  the  previous  expression  for  the  conditional  probability 
density  functions  which  are  derived  assuming  a  'mark'  signal 
has  been  sent.  That  is,  for  a  given  value  of  q-^,  an  error  is 
made  if  qg>q-^.  Thus  the  average  error  probability  is  found 
by  averaging  the  conditional  error  probability  given  by 

Pel  =  Pr<  *€>  >  <*1  I  Hl  > 
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GO  i*     GO 


p       IU  p       HI 

"     [     <q*'V  dq®  I  PQ1<ql|Hl>  dql 
J  0   Ja. 


(3.12) 


Substituting   Equation   (3.11)   in   Equation   (3.12)   and 
interchanging  the  order  of  integration,  we  obtain 


J0   J0    Jc 


q0  *m01 


l0      2  o 


qoy  m0i       1 

•Io(     2    )>dq0    dff 

a  ' 

PQ  <q1|H1)  ]  dqx     (3.13) 


In  the   above   equation,   the   inner-most   integral   can  be 
expressed  in  terms  of  the  Marcum  Q  function  as  follows 


f 


q0  +In01 


<\o/~™r 


Jrnl 


l0      2   a~  ,H0"  m0l  ,  ^TO0i   ^i 

•I0( j )  dq0  =  Qc  -|i.-J-  > 


Then,   Equation  (3.13)   for  Pen   with  the  aid  of  Equation 
(3 . 10 )  becomes 

D»2n         j 

-n-r 


m0i    qi 

Q(    — . )    d* 


■  f-4-  r  a- 


a    '    a 


ql+mll 


2   a2      _    #V^TI 
V 2— 


)    d0 


]dql 


(3.14) 


From  the  orthogonality  property  of  the  signal  pair  used 
(which  is  obtained  by  assuming  sufficient  separation  between 
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the  two  frequencies  and  that  c^i    as  well  as  coq    are  large), 
we  have 


f 


sin<&>  +t  )  •  sin*>  t  •  dt  =  0 


so  that  the  term  itiq^  is  independent  of  0      Therefore  Equation 
(3.14)  can  be  rewritten  in  the  following  form 


el 


h  "i  r 


'Sal   qi 

9<  — — .— -  ) 


qi+mn 

2    a2 


qi,/ivr 

Io<   J— T1  >   dq 


d£ 


(3.15) 


where  the  order  of  integration  has  been  changed. 
Furthermore,  using  the  following  formula  involving  an 
integral  of  a  Marcum  Q  function  [Ref  5] , 


r 

Ja 


a.    *V 


2   2   al 


«  V 

V  -S  )dV 


2X   2 
al  +<T2 


1  -  Q 


2,   2' 
al  +a2 


2A   2 
al  >a2 


<7. 


2x   2 
al  +a2 


»  a.    +a_   >l   a.  +<y 


the  inner  integral  in  Equation  (3.15)  can  be  simplified  in 
such  a  way  that  P  ^  becomes 

r2n    r- 


.-i 


1  - 


2n 


Q< 


^^7 


)  d0 


<d 


2  a  V2   a 
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(3. 16) 


where  the  dependence  on  0    is  imbedded  in  the  term  m-^  only. 

Following  exactly  the  same  procedure  used  in 
obtaining  the  expression  for  Pei>  it  can  De  established  that 
the  expression  for  P  q  which  denotes  the  error  probability 
when  Hq  is  assumed  to  be  the  true  hypothesis,  that  is, 
Pr(q-^>qQ  |  Hq  )  takes  the  form 


,  if  i  -  _»_r 

»0      2  l  2n 


_2°_,_i°_>  d. 


i  r2"    ^ 


y; 


to 


00 


wT 


)    d* 


<3.17> 


where  only  the  term  hiqq  is  dependent  on  6  .  Therefore  the 
total  average  Probability  of  error  Pe  can  be  obtained  from 


■  Hh  |  Pel  *  %»  J 


assuming,  as  previously  indicated,  that  the  two  hypotheses 
are  equally  likely.  Using  Equation  (3.16)  and  Equation 
(3.17),  it  is  obtained  that 


-  <l 


d0 


QC-22-,- 


i°->  -  Q(-ig-,   gg  ) 


H 


(3. 18) 
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If  we  now  use  the  jamming  waveform  which  is  optimum 
under  energy  constraint,  when  coherent  FSK  reception  taken 
place  that  is 


n.(t)  ■ 
J 


nj 


4      T 


Csin«»ft  -  sintf.tl 
l         io 


(3.19) 


then  the  terms  m^^(i,j=0,l)  in  Equation  (3.18)  which  are  a 
function  of  the  jammer  waveform  n^(t),  can  be  computed  as 
follows 


™11  s  C(ASi,«,Sl>  +  (nj,Sin2  *  C<flSii5*Ci)  *  <n.^C,)22 


flpp      2            ATv1r"~ 
=    (   flI   )      ♦  !li. 


cos9    + 


P     .T 
nj 


(3.20) 


ro_       a    C ( AS 
01 


P     .T 


i,i'V  *  <v  V1    +  c<ASi,**V  ♦  <nrc0n 


(3.21) 


and 


m00  =   C(AS0,a'V    +    <VV3      +   t(AS0,*'V    ♦    c«rt#u: 
!  ATSF 


=  <-SX->    . 


nj 


cosO    + 


P      T 


(3.22) 


mlO  "  C(AS0f5»Si)  +  (VS1)3   +  L<ASV,B'C1>    *    <»..Ci>J 


P  .T 
-   nJ 


(3.23) 


Thus  the  Probability  of  error  can  be  expressed  in  terms  of 
SNR  and  JSR  only,  as  follows 
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p2n 
J© 


11 


'a 


01 


^2     /t 


)    -    Q< 


01 


10 


)    -   Q( 


10 


00 


yr     yr 


dtf 


(3.24) 


where 


axl  =  SNR-(2  ♦  2VJSR-cos0  +JSR) 
*01  =  al0  "  SNR'JSR 


a__  =  SNR-(2  -  2,/tfSR-cos*  +JSR) 


SNR  = 


JSR 


2 
AT 

2 

No 

Receiver  performance  can  now  be  evaluated  as  a  function  of 
SNR  for  fixed  values  of  JSR,  using  Equation  (3.24). 

These  results  have  been  included  here  so  that  it 
will  be  possible  to  compare  receiver  performance  given  by 
Equation  (3.24)  with  performance  of  a  similar  receiver  under 
the  assumptions  of  jammer  miss-synchronization  or  frequency 
offsets.  In  every  case,  it  has  been  found  that  a  degraded 
jammer  produces  a  smaller  receiver  Probability  of  error  than 
a  'perfect'  jammer.  This  is  to  be  expected  and  the  analysis 
of  the  next  two  sections  demonstrates  this. 
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C.   ANALYSIS  OF  THE  EFFECT  OF  THE  JAMMING  WAVEFORM  HAVING 
RANDOM  TIME  OF  ARRIVAL 

1.   ASK 

In  the  previous  section,  the  effect  of  the 
deterministic  jamming  waveform  given  by  Equation  (3.8)  on 
the  performance  of  the  ASK  receiver  of  Figure  (5.3)  was 
studied.  Since  such  an  ideal  jammer  would  have  to  be 
synchronized  with  the  signal  transmission,  we  now  develope  a 
more  realistic  model  in  that  the  difference  in  the  arrival 
time  of  the  jammer  waveform  with  respect  to  the  signal  can 
be  characterized  by  a  parameter  t  •  Thus  the  difference  in 
the  time  of  arrival  between  the  jammer  waveform  and  the 
information  waveform  is  given  by  the  parameter  r  so  that  we 
have 


n.(t)  . 

J    t 


2 
nj  T       c 


The  effect  of  having  this  time  difference  t  will  now  be 
analyzed  using  the  results  of  the  previous  section. 

First,  using  Equation  (3.1),  for  a  given  value  of  t 
the  conditional  means  of  X  and  Y  become 


E<X|Hlt*>T  =  (flS^.S)  ♦  <n.,S>T  =  ^|Hii,iT 
KY|Hlf«T  =  <ASe.C)  ♦  <„..C)t  =  -,,„   ,(T 


Since  the  conditional  variance  and  covariance  of  X  and  Y  are 
independent  of  the  signal  and  jammer  waveform,  the  results 
on  the  variance  and  covariance  of  X  and  Y  derived  previously 
in  section  B.l  can  be  used  here  to  obtain,  for  any  given 
value  of  t 


Var<X|H  ,*>   =  VarCY|Hif0>   = 
xt  'It 


NoT 
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The  mathematical  expression  previously  derived  for  Pe  is 
useful  here  in  order  to  obtain  the  ASK  receiver  performance. 
Some  modifications  however  are  necessary  that  require  the 
computation  of 

(ni»S)*-  =    JP«;-JF-  'Sinn  (t-T>«sin»  t«dt 


,s)  =  [ft? 

j'   t      4    nj  T 
J© 

T r~      r  2«sin«#  T-cos(»  T-»  t)   , 

■  JPnj4-    ««V  *  ~ °    ° 

1  2sin«  T  J 


=  Jp„ i — o —  'COStt  T 

v    nj  2         c 


Since  wcT=n<jr   ,   the  second  term  in  the  above  bracket  is 
zero.   Similarly,  we  can  show  that 


(n.,C)   =   P  .-^—  •sin(-«  t) 
j    t    4    nj  2  c 


These  results  are  used  in  the  evaluation  of  Equation  (3.2), 
namely 

.2 „,   .  .2^ 


=  C(AS^,S)  +  (n^S)^.]2  +  C(AS^,C)  ♦  <njfC>Tl 

=  <***  >   +  at|p   4-'co5(V+5)  *  P«i-T-        <3.25) 
2         m  nj  2        c         nj  z 


and  in  the  evaluation  of  Equation  (3.3),  to  yield 

M»   =  (n  ,S)T2  ♦  lnycyr2 

T  (3.26) 

=  P  -=r- 
nj  2 

Using  Equation  (3.7)  along  with  the  results  of  Equation 
(3.25)  and  (3.26),  for  a  given  value  of  j  it  can  be  seen 
that 
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P  (T)  = 

e 


r2"  JrT  Jfv"        i 

*  J  A 


(3.27) 


which  in  turn  can  be  expressed  in  terms  of  SNR  and  JSR 
Since 


M 


-J—  -  2SNR  ♦  4SNR«^JSR-C0S<«  t+0)  ♦  2SNR-JSR 
2  o 


and 


M1 


»  2 -SNR -JSR 


substituting  the  above  in  Equation  (3.27),  P„(t)  becomes 


P  <T>   = 

e 


i  ♦  Q(,/2SNR-JSR,-li-> 

1  o 


I 
2n 


p2n 

Q(J2SNR(1  +  2,/TSRcob(6»  r*6  )  +JSR  )  .  -^>  .  dff  1   (3.28) 
0 


Observe  now  that  the  integration  of  Equation  (3.28) 
involving  the  variable  0  is  over  the  range  [0,277]-  Thus 
regardless  of  the  value  taken  by  the  term  o>c7  >  since 
cosine  is  periodic,  the  values  taken  on  by  cos(cj  7  +  0)  are 
the  same  as  the  values  taken  on  by  cos  Q  ,  with  O<0<2tt 
Thus,  Equation  (3.28)  is  in  essence  identical  to  Equation 
(3.9)  which  yields  Pe  when  the  jammer  waveform  has  no  timing 
difference  7  .  Therefore  it  is  clear  that  any  timing  error 
associated  with  the  jammer  waveform  doesn't  affect  the 
performance  of  the  incoherent  receivers  for  ASK,  when  the 
jammer  waveform  is  given  by  Equation  (3.8). 
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2.   BFSK 


In  this  section,  an  analysis  similar  to  the  one 
carried  out  in  the  previous  section  is  now  applied  to  FSK 
modulation.  The  timing  error  associated  with  the  jammer 
waveform  for  the  FSK  system  is  modeled  by  modifying  Equation 
(3.19)  as  follows 


nj(t)r  "  JP„j4-  •[•i-Vt- 


t)  -  sin«_(t-T)3 
0 


(3.29) 


As  in  the  previous  section,  it  is  found  that  changing  the 
jammer  model  requires  only  a  recomputation  of  the  terms  m-t-t  , 
m01»  m00'  anc*  m10  £iven  by  Equations  (3.20),  (3.21),  (3.22) 
and  (3.23)  respectively.  Thus,  using  Equation  (3.29),  we 
will  compute 


NoT 


It  can  easily  be  seen  that 

VarCX1|H1,fl>T  «  Var<YJ[  |H1 , 0>t  = 

so  that  these  variances  remain  unchanged.  Furthermore, 


n.,S.)   =    |p   -i 
j*  1  t      4    nj  T 


Csin«t(t-r)  -  sin»_ ( t-T ) 3 • sin«, t «dt 

1  0  1 


j  2   -cos«iT 


and 


(n.fC.  ) 


Jp  .-4-   tsin«,(t-T) 

<  nj  T         1 

> 
J  nj  2 


sin«.( t-T) 3  cosw.t «dt 
0  1 


sin(  "**,▼ > 
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Therefore,  from  Equation  (3.20),  we  obtain 


'ilT  "  C<ASli^»Sl)  +  (nj'Si>T1    +  C<ASl,«'Ci)   *  Cn4'Cl}T] 

2     ATVP" 


<J^>    ♦ 


P   T 

^= COS(«*1T+5)   ♦  ^— 


and  from  Equation  (3.21),  we  obtain 


P  .T 
nj 


01 T 


Also  from  Equation  (3.22),  we  have 

J    0  T 


«  C(AS0  ^,S0)  ♦  C«irB.>_l?  ♦  C<AS„  fl,Cft)  ♦  <n  .  ,Cft>„:i2 


00T 


0,0*  0 


V^0  T 


._   2     AT^P 


P  .T 


and  from  Equation  (3.23),  we  have 


m 


i*r  =  c<AS0f^si)   +  (nj'si)r32  +  C(flS0,a»ci)   *   (nj'ci}T32 


P   T 
nj 


By  substituting  the  above  results  in  Equation  (3.24)  we 
obtain  Pe  for  a  given  value  of  t  namely 


'A « -^ 


r2n 


0 


q<  _LL_,_°JL  )  -  q<  —^L-.—Li 


/T 


/T         /~2 


4n 


r2n 


0 


y~2         S~2  S~2         S~2 


]•**] 


where 


a         =   SNR«(2    +    2v^SR-cos(«*1T+0)    +JSR) 
«01    =    *10   "    SNRJSR 
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oM  m   SNR-<2  -  2,/JSR-cos<»-T+#>  +JSR) 


SNR  = 


2 
AT 

2 

No 

P 
nj 

JSR  = 


By  observing  above  probability  error  equation  it  is 
noted  the  the  only  term  changed  is  cos#  .  Although  the  cosfl 
term  is  changed,  due  to  the  integration  of  the  Q  over 
interval  [0,2tt],  whatever  the  value  j  is,  the  Pe  doesn't 
change.  The  timing  error  associated  with  the  jammer  waveform 
does  not  affect  the  performance  of  the  incoherent  receiver. 

An  observation  similar  to  the  one  made  at  the  end  of 
the  previous  section  reveals  here  also  that  Pe(r)  is 
actually  independent  of  t  •  Therefore,  timing  errors 
associated  with  the  jammer  of  Equation  (3.19)  do  not  affect 
in  anyway  the  performance  of  the  receiver  analyzed  in  this 
section. 

D.   ANALYSIS  ON  THE  EFFECT  OF  THE  JAMMING  WAVEFORM  HAVING 
RANDOM  FREQUENCY  ERRORS 

In   this   section,   the   effect   of   frequency   errors 

associated  with  the  jamming  waveform  is  analyzed  insofar  as 

the  performance  of  incoherent  binary  receivers  is  concerned. 

Results  derived  in  section  B  that  are  applicable  to  the 

present  problem  will  be  utilized  in  this  section.  As  is  will 

be   seen,   frequency   errors   associated   with   the   jamming 

waveform  change  only  the  mean  value  of  the  receiver  output. 

This  fact  makes  the  analysis  of  the  receiver  performance 

more  tractable. 

1.   ASK 

In  the  previous  section,  we  studied  the  effect  of 
timing  errors  in  the  jamming  waveform.  Now  we  assume  that 
the  frequency  error  of  the  jammer  waveform  is  Ao>  while  the 
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timing  error  is  zero  or  can  be  shown  to  be  negligible.  Thus 
we  have 


nj(t)    =  K  sin(o)c  +  Aco)t 

Since   the   jammer  power   is   constrained   to   be   pnj  > 
constant  K  must  be  consistently  defined.  Since  we  have 


the 


..2      -2     2, 
n  .11   =    K  >sin  (« 

J 

Ja 


+A«)t  dt 


2 

K  T 


we  must  have 

K  ■ 
where 


1  - 


sm2(i»  +  d**>T 
c 

2(&v  +  /to)T 

c 


nj 


4      T-CK»] 


K 


1  - 


sin2(6»  +A^)T 
c 

2(w  +  4«>T 
c 


=  P 


nj 


In  order  to  evaluate  the  receiver  P  using  Equation  (3.1), 
for  a  given  value  of  Acj  »  the  conditional  means  of  X  and  Y 
become 


E<X|H  ,0> 


/W 


(AS_,S)  ♦  (n.,S) 


j'   44*    mK\Hlt6,Au 

1  1 "   4«        0 '         J '   A«     Y  J  H  ,  8 ,  At* 

where  Sn  C,  and  S  have  been  defined  in  the  previous 
chapter.  It  can  be  shown  that  the  variance  and  covariance 
of  X  and  Y  conditioned  on  Hi  and  Q  are  independent  of  the 
signal  and  jammer  waveform,  so  that  in  spite  of  the 
frequency  error  Aio  present,  we  have 


VartXlH. ,8) A      =  VartY |H. , 0} 


NoT 


Au 


and 
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EC  CX  - 


■klH^.Ai1;"  -  "YlH^O,**1     l"i^^  5  .  O 


In  order  to  compute  the  mean  value  of  the  receiver  output 
under  the  assumption  of  jammer  frequency  error  as  described 
above,  we  must  compute 
rT 


j'  Au  c  c 


r 

Ja 

sinAt>T  _  c 

d<*T        (2m  +iW>T  ] 


0 

sin(2w  +A«)T 
KT_ 

^»T        (2« 
c 


and 


<n.,C).   =    K-sin<«*  +A»)t'Cos«*  t-dt 


r 

■    K*sin 
J® 

1  -cosAaT  c j 

AuT  <2«  +/W)T  j 


1  -cos(2<*  +A<*)T 
KT 


With  these  results,  by  substituting  in  Equations  (3.2)  and 
(3.3),  we  have  the  conditional  means  M  and  M',  given  by  Ma, 

and  m'a« 

AT 
H,  =  (-^->  ♦  ATECn-.S).  'cos©  +  (n..C)»sin0] 

Au  2  j    Aa  j ' 

♦  (n.,S)  .  2  ♦  (n.,C).  2 
j    At*       j    A« 

AT  2      AT  9   9       9 

=  (-4p->   ♦  _2l_.K.TLB»3  +  4K  -T   CC»3  (3.30) 


where 


[sin(A«T-fl)+si 
=  [     A«*T 

I   A«*T((2»  ♦A«)T>  ' 


sin((2«*  +A«)T  +  sintf 
n#  c 


(2<*  +A<*)T 
c 


and 
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M»    «  (n.,S)  .  2  ♦  (n  ,0.  2 

At*  J '  Au  j '  Au 


^.T2- 


EC'l 


(3.31) 


By  substituting  the  above  results  in  Equation  (3.7),  for  a 
given  frequency  error  Aco  5  the  receiver  error  probability 
becomes 


P  <4«> 


Jo 


*^-,-£)d0  +Q( *• 


■•*! 


(3.32) 


Equation  (3.32)  must  to  be  expressed  in  terms  of  SNR  and 
JSR.   Since  it  can  be  shown  that 


M 


Ac* 


a 


2 

A(A 


=  2SHR  +  4SNR->/jl5R- 


CB'3 


EK»3 


1/2 


♦  32 -SNR -JSR 


EC*  3 
EK»3 


=  32  SNR -JSR 


■  T2 


EC'3 
EK'3 


where 


SNR  =  AzT/No 


JSR  =  Pnj/E 


Equation  (3.32)  now  becomes 


P    (At*)    = 


1    +    Q( 


32  •SNR- JSR       CC,]     ,-5— ) 
EK'3  a 


r2ir 


1 

2n 


Q< 


a 


2SNR+  4SHJU/JSR EB'  I  i„+32SNRJSRCC  *  3     .  -2-)  dtf 


EK»3 


1/2 


EK«3 


(3.33) 


where  the  dependence  on  A  CO  is  imbedded  in  the  terms  B',C 
and  K ' . 
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By  defining  a  normalized  frequency  error,  Acjjj  >  the 
above  results  can  be  rewritten  as  a  function  of  Aco^»  where 

N     a 


The    terms    B',C'    and   K'    as    a    function   of   A^n    Decome 


B» 


C» 


K» 


■1 


sin(tt   T* A«„-0>+sin0 
c N 

c  N 


Am 

s  i  na   T • - 
c  l 


N 


.  f !_^ 1 

Ac*.  ■«    T(2+^6>-,) 
1  N       c  N 


sin26>   T(l+ite„) 
c N_ 

c  N 


sin(t»   T(2+A<a„>    +    sin0 
c N 

«    T(2+A«„) 
c  N 


so    that 


e         N  2 


1    ♦    Q( 


32SKR-JSR 


[C'3 
[E»  J 


r2n 


2n 


Q< 


0 


2SNR+4SNR^JSR [B'?  _  +  32SNRJSR-C ( 


CK»  ] 


1/2 


CK»3 


>d0 


(3.34) 


The  developed  expression  is  actually  the  probability 
of  error  conditioned  on  a  given  frequency  error.  A  procedure 
similar  to  the  one  used  to  analyze  the  effect  of  frequency 
errors  on  a  coherent  receiver  is  used  here  in  order  to 
obtain  the  average  probability  of  error  for  the  incoherent 
receiver  of  Figure  (5.3).  It  is  clear  that  the  conditional 
probability  of  error  for  a  given  value  of  Ao>m  is  a  function 
of  AcOfj.  Therefore,  in  order  to  compute  the  average  value  of 
Pe ,  it  is  assumed  that  Acom  can  be  modeled  as  a  random 
variable  whose  probability  density  function  is  uniform  over 
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some  range  (AconL'^^NH^  »  where  Ao>NL 


is  some  the  lower  limit 
and  AcdpjH  is  some  upper  limit  of  the  normalized  frequency 
error.  With  this  assumption,  the  average  Pe  can  be  obtained 
from 


P   = 

e    A& 


Am 

I  r   M,H 

M,L 


P  (Aw.,)  'dAu.. 
e    N      N 


or  equivalently 


At* 


1  +  Q< 


32-SNR-JSR  JL£11  ,-£_> 
CK»J    ° 


j  p      N,H 

H,L 
an      -. — 

Q(    2SNR(1*2,/7SR CB'  ]  /0  +  16  ■  JSRCC  '  3     ,JI-)dg]    dAtt„ 

«    J  CK'31/2        CK'3    a     i     N 


_1_ 
2n 


9 


(3.35) 


where  B',C'  and  K'  have  been  previously  defined.  The 
computation  of  P  is  complicated  by  the  fact  that  B',C'  and 
K'  are  themselves  function  of  Aco^' 

2.   BFSK 

In  the  previous  sections,  a  deterministic  waveform 
has  been  used  to  model  a  jammer  signal.  Here,  we  introduce 
a  more  realistic  jammer  model  which  includes  a  frequency 
error.  The  reason  for  choosing  such  a  model  stems  from  the 
fact  that  in  practical  situations,  a  jammer  will  lack  exact 
knowledge  of  the  transmitting  frequencies.  Thus,  for  BFSK, 
we  model  the  jammer  as 

j    Au    ,A<*  11  0    0 

Due  to  the  power  constraint  on  the  jamming  signal,  we  must 
satisfy 
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Hn.<t>H    2  =      K2    fsin<«    +  4«,  >t    -    sin(«_-n4»_ )  J2  -dt 


■■"•"J, 

i- 


_.2_       |              sin2(«t+Attt  )T              sin2  <**    +/W_  )T 
S    T         |  „  1 1 ^    0 0 

2{<*.+A*.)T  2(«_+4«_>T 

11  19  V 


sin<«».+Aa .  )T  sin(«    +  iW    )T       , 

-       2 2 ^ ♦       2 ^ 5 =      P        : 

l».*Au.)T  <«    +A«    )T       '  nJ 

d        d  s         s 

This  means  that  the  jammer  waveform  must  be  reexpressed  as 


IP         2 


n  (t).        «   -s Csin(»t+A»t  )t  -  sin(«_+A«)tl 

j    Au    ,A(*9  1  11  « 


2 
CE-3 


where 


K"  ■    2-  sinc2(»1+A«1>T  -  sinc2 laQ+Aa^)T 


-    2sino(».+ii»  -  >T  +  2aino(«#  +A«*  >T  (3.36) 

ad  s    s 


and 

sin  X 


X 

s              10 

A«  =  A&.+A&- 
s               10 

d            10 

4«  a  At* .-An- 
ti              10 

By  assuming  that  the  jammer  waveform  contains  frequency 
errors,  it  turns  out  that  the  only  changed  statistical 
parameters  are  the  mean  values  of  the  receiver  output 
conditioned  on  a  given  frequency  error.  Thus  we  can  use 
Equation  (3.18)  in  order  to  compute  Pe  as  a  function  of  the 
frequency  error.  We  need  to  compute 

rT  _J_ 


w*v4.0  -j  [  [  VZ 


(sin(M1^^£»1)t  -  sin<«  ♦*«  )t>  I -sim*  t  ldt 
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=    ]p         -JL_  -CK-3       2CB-3  (3.37) 

>l    nj      T 


where 


sinA«.T      sin(2», +il».  )T      sin(».-A»_)T      sin(«    +4«_)T 
B~    _   1_  _ 1  1 d         0         + s         0 

Att.T  (2c*f+A<«,)T  («*    -«_)T  (4»    +ii«_)T 

1  11  d      0  s         0 


(3.38) 


Similarly,    we    can   show   that 


ji^         >l      nj      2 


aj»Ci>^ 


where 


l-cos(2«1+A«1)T  l-cos4»   T  l-co5(#s+4«0)T 


C"    ■ + 


(2«    +il»1)T  ^T  (»s*i4»0)T 

l-cos(«  .+4«_)T 
a        w 


(3.40) 


(«.+4tt_)T 
a         © 


Also 


"    »s«>a         a         «    Jp„4    -i-CK-l-tD-l 
j'    0   Au.,Aa_         >l    nj      2 


(3.41) 


where 


*in4»0T      sin(2«0-*46»0)T      sin<<*    ♦£«    )T      sin(«    +4»    )T 
A<*qT  (2«0*iSc*0)T  («d+i4»1)T  («s+/l«    )T 

(3.42) 


Finally 
( 


Wa^.a.    ■  J  pnj  4-  ;,c*"]    2cfm:i  (3'43> 
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where 


l-cos(«  ,+Aa. )T 
1 +  5 1 

<U.+At*.   )T 

a    l 


1 

-cos  (a    +At*.  )T 
s          1 

♦ 

(«*    +Aa.  )T 
s         1 

♦   . 

co  s ( 2»_ + Aw _ ) T 
so        0 

cosAu.T  -  1 
♦  (3.44) 


(2<V*VT  *v 


Using  Equations  (3 . 20 ) , (3 . 21) , (3 . 22)  and  (3.23),  the 
conditional  means  of  either  channel  output  conditioned  on 
either  hypothesis  and  conditioned  on  a  given  frequency  error 
can  now  be  computed.  From  Equation  (3.20),  we  have 

m, .        =  [(AS,  ofS.)  +  (n.,S.),    .  I2 
11.    .  1,5*1       j '  1  Aa ,  ,  Au-. 

A<* .  ,  At*  '  w        1 '   0 

♦  t(AS,  .,C,)  ♦  (n.,C. ) .     .   32 
1,8'     1        j'  1  A»lfA«0 


2 


-?^)  cos2©  +  2(-^I)cosd.  P  .  -£ 
2  2       ^  nj   2 


2 
T      -1      2     AT       2 
+P    -JL.£K~J    CB"3   +  (=^— )  sin  8 
nj   2  i. 


=  <£*_>   + 


AT  T  2  T         - 1 

-2<-^4->sin0  P    4—  CK"1    CC"]  +  P    -J — CK"3    CC"3 

2      ^  nj  2  nj   2 

1^ 

'J?  .J— CK"]  2.f  CB"]cos0  +  CC-]sin0 

_1_ 

+  P  .-?-•  CK"3  2.  I  CB"32  +  [C"]2 
nj  2  I  J 

From  Equation  (3.21),  we  have 

m_,        =  [(AS.  -,S_>  +  (n.,S_).       32 
01.     .         1,0*  0       J   0  Au..Au- 
At*    ,Au  '  ■        1'   0 

♦  [(AS.  „,C_)  ♦  (n.,C-),       I2 
1,00       j  '  0  At*  t  ,  ^«„ 

1      17 

T        -1       2         T        -1       2 
=  P   -J—  [R-3   .[D"3^  +  P   _i_.[K"l    [F"l 
nj  2  nj  2 

Also  from  Equation  (3. 22), we  have 
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At*    ,Au  '  ■        1*0 

♦  C<AS_  -,C_)  «•  (n.,C_).    .   32 

0,0*  0       J   0  ^.i^. 

2       1 

-  <^->   ♦  AtIp  Z—  CK"3  2-|  CD"]cos0+CF"]sin0 

2  \i  nj2 

1 

T        2  f      2      2 
♦  P  .-i—'£K~3   •   CD-I  +CFM3 
nj  2 

and  from  EQuation  (3.23),  we  have 

m  -  C(AS_  -,S.)  +  (n.,St).    .  I2 

♦  CAS_  -.C.)  ♦  (n^C.)  .    .  I2 

0 1 S '    1       j   1  /W ,  ,  4u  _ 

1     %9 

T        -1       2         T        -1       9 
-  P   -75 — CKM1   -IB"!'  ♦  P  .-i—IK"]   .£C"1 
nj  2  nj  2 

With  SNR  and  JSR  as  is  previously  defined,  we  can  rewrite 
the  above  equations  specifying  the  conditional  means,  as 

1 

rai  *  r  t o 

( ii)  A  m    2SNR-       1    +2v/JSRCK"3     (CB-3 -cos^+CC-] -sin») 

a2   A"i'^0        I 

♦  JSR-tK-]~1(CB,,32+CC"32) 

( — ®4>  -       =  2SMR-JSR-tK"3"1' (CD"]2^CF-J2) 
a     i*   0 

-i 
( — ?®->  .        »  2SNR-  I  1  ♦2^JSR-CK"]2(CDM3 -cos^  +  CF-3 -sin») 


a2   *V**0 


♦  JSR-CK-]  1(CD-]2+tF"a2) 


( ~)  .       =  2SNR-JSR.tK"3_1- (CB"32+CC-]2) 

a2  ^,^0 


Using  the  above  results,  we  obtain  the  conditional 
Probability  of  error  by  substituting  in  Equation  (3.18)  to 
yield 
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p  <4«,»4«*>= 


e if   0    2 


where 


X 

alx  =  2SMR- I  1  +2^JSR-CK~32(CB"1 -cos^+CC"! -sin0> 
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+  JSRCK"j"1(CB"32*CC"32> 


(3.45) 


1  (3. 


46) 


a^,  =  2SNR-JSR-CK-3"1- (CD-32+CFM32) 
vl 


(3.47) 


a.„  =  2SMR.JSRCK"!""1- (CB"]2*CC-]2) 


(3.49) 


As  in  the  previous  chapter,  normalized  frequency 
errors  Acom  i  ,  and  Acom  q  are  introduced  and  it  is  assumed 
that  these  normalized  frequency  errors  can  be  modeled  as 
independent,  uniformly  distributed  random  variables  over  a 
given  range.  Thus,  following  the  same  procedure  previously 
used,  the  average  probability  of  error  can  be  obtained  as  a 
function  of  Acojj  i  and  Aco^  q*  Using  the  assumption  of 
independence  and  of  uniform  distribution  of  the  frequency 
errors,  we  obtain 
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e        A<* 


1  1  r*iHillHrA"lll0fII 


^"n.I.L   ^"n,©,!- 


V*"N,l»*"Nf0)dilaN,ld^N,0 


Using  now  Equation  (3.45),  we  obtain 


e~   A« 


H, l,H"4wNf 1,L 


^•••H"^Mf0,L  J         L 


J® 


U0 


d*WH,ldi4wN,0 


(3.50) 


where  a-^,  q^  >  «3 »   anc^   «4  have  been  defined  by  Equations 

(3.46) , (3.47) , (3.48) ,   and   (3.49)   respectively,  and   the 

parameters  K" ,B",C"  and  F"  have  been  defined  by  Equations 
(3.36), (3.38), (3.42)  and  (3.44),  respectively. 
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IV.  DESCRIPTION  OF  GRAPHICAL  RESULTS 

A.  GENERAL 

In  previous  chapters,  mathematical  models  of  optimum 
jammers  were  introduced  in  such  a  way  as  to  obtain  more 
realistic  jammer  models  which  include  timing  errors  and 
frequency  errors.  Analyses  were  carried  out  in  order  to 
determine  the  effect  of  the  timing  errors  or  frequency 
errors  on  the  performance  of  coherent  as  well  as  incoherent 
digital  receivers.  This  work  has  resulted  in  four  main 
equations  which  specify  the  performance  of  the  receiver  in 
terms  of  the  probability  of  error  as  a  function  of  signal  to 
noise  ratio  (SNR)  and  jammer  to  signal  ratio  ( JSR) .  With  the 
aid  of  computer  evaluations,  the  receiver  probability  of 
error  for  various  jammer  models  has  been  obtained  as  a 
function  of  SNR  for  a  given  JSR  value  and  the  results 
plotted  for  each  case  considered. 

In  each  plot,  the  JSR=0  case  has  been  included  in  order 
to  provide  a  basis  for  comparisons  of  the  jammer 
effectiveness  on  the  receiver  performance  as  it  relates  to 
additive  white  Gaussian  noise  only  interference. 
Additionally  the  effect  of  a  perfect  jammer (t-0  ,Aco  =  0)  is 
also  evaluted  in  order  to  determine  jammer  degradation 
resulting  from  timing  or  frequency  errors. 

The  effect  of  the  jammer  waveform  on  the  coherent 
receiver  will  be  presented  first  and  then,  the  results  for 
the  incoherent  receiver  will  be  presented  next. 

B.  COHERENT  RECEIVERS 
1.   PSK 

Graphical  results  on  the  performance  of  a  coherent 
BPSK  receiver  in  the  presence  of  optimum  jamming  are 
presented  in  Figure  (5.6),  based  on  the  results  given  by 
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Equation  (2.7).  Plots  of  Pe  were  computed  as  a  function  of 
SNR  for  fixed  values  of  JSR. 

Figure  (5.6)  clearly  shows  the  'break  point' 
phenomenon  described  in  [Ref  1]  ,  namely  that  as  a  JSR 
increases  to  a  value  of  one  or  greater,  Pe  does  not  decrease 
with  increasing  SNR.  That  is,  for  JSR  value  greater  than  or 
equal  to  one,  Pe  increases  to  the  value  of  1/2  with 
increasing  SNR.  From  this  figure,  it  can  be  noted  that  10.2 
dB  of  SNR  is  required  to  obtain  a  Pe  of  10" 6  at  a  JSR  value 
of  0.  In  comparison,  it  takes  14  dB  of  SNR  to  obtain  the 
same  Pe  for  a  JSR  value  of  0.1. 

Figure  (5.7)  displays  the  probability  of  error  for 
the  same  coherent  receiver  under  the  assumption  that  the 
optimum  jammer  previously  considered  suffers  from  a  timing 
error  xw  which  is  assumed  to  be  uniformly  distributed  over 
the  inteval  [0,1]  . 

It  can  be  seen  that  now  13  dB  of  SNR  is  required  to 
obtain  the  Pe  of  10  °  at  JSR  =  0.1.  We  can  compare  this 
result  to  the  14  dB  SNR  required  when  the  jammer  model  does 
not  include  a  timing  error.  It  is  apparent  however  from 
these  results  that  the  penalty  for  lack  of  coherence  on  the 
part  of  the  jammer  with  the  signal  transmission,  is  not 
severe.  Therefore,  without  knowledge  of  correct  timing  for 
synchronization  with  the  adverary's  coherent  receiver,  the 
jammer  has  to  increase  its  power  by  a  small  amount  in  order 
to  produce  the  same  error  rate  at  the  receiver  when  jammer 
synchronization  is  perfect. 

When  the  jammer  model  incorporates  frequency 
estimation  errors,  the  performance  of  the  receiver  is  given 
by  Equation  (2.33)  and  Figure  (5.8)  shows  the  resulting  P 
of  the  receiver.  The  plots  shown  in  this  figure  were 
obtained  under  the  assumption  that  the  normalized  frequency 
error  A60jj  can  be  modeled  as  a  random  variable  uniformly 
distributed  over  the  range  which  in  this  case  has  been 
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assumed  to  be  [-0.001,0.001]  and  that  the  number  of  waveform 
cycles  in  one  observation  period  is  10.  By  analyzing  Figure 
(5.8),  it  is  noted  that  Pe  of  10"6  can  be  obtained  with  13.6 
dB  of  SNR,  when  JSR  is  0.1. 


TABLE  1 

SNR  REQUIREMENTS  (P£  =  10 "6  PSK) 

JSR 

DET'  JAMM 

JAMM  WITH 
TIME  ERR 

JAMM  WITH 
FREQ'  ERR 

JSR=0.0 

10.2  dB 

10.2  dB 

10.2  dB 

JSR=0.1 

14.0  dB 

13.0  dB 

13.6  dB 

JSR=0.5 

26.0  dB 

20.0  dB 

20.0  dB 

JSR=1.0 

unable 

large 

46.0  dB 

The  values  in  Table  1  were  obtained  directly  from 
the  figures  corresponding  to  Pe  when  the  jammer  exhibits 
timing  errors  or  frequency  errors. 

2.   BFSK 

The  procedure  used  in  the  previous  section  is  also 
appropriate  in  the  analysis  of  the  figure  presenting  the 
results  on  the  performance  of  the  coherent  BFSK  receiver. 
Figure  (5.9)  displays  the  P  of  this  receiver  in  the 
presence  of  an  optimum  deterministic  jammer  which  exhibits 

no  timing  or  frequency  errors.  One  can  note  that  when  JSR  = 

-  6 
0,  13.4  dB  SNR  is  required  in  order  to  obtain  a  Pfi  of  10   , 

whereas  when  the  JSR  value  increases  to  0.1,  18  dB  of  SNR  is 

required  to  maintain  a  Pg  of  10"  .   Figure  (5.10)  shows 

plots  of  Pe  for  the  same  coherent  receiver  when  the  jamming 

-  6 

waveform  suffers  from  timing  errors.  For  Pe  of  10   '  at  JSR 

=  0.1,  17  dB  SNR  is  required.  This  is  1  dB  less  SNR  than 
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that  required  for  the  case  of  a  jammer  waveform  with  no 
timing  errors.  In  order  to  obtain  these  results,  it  was 
assumed  that  the  normalized  timing  error  t^>  can  be  modeled 
as  a  uniformly  distributed  random  variable  over  [0,1]. 

Figure  (5.11)  displays  plots  of  Pe  for  the  same 
receiver  when  the  jammer  is  assumed  to  contain  frequency 
errors.  This  figure  was  obtained  under  the  assumption  that 
n=6,  m=4  where  'n'  is  the  number  of  cycles  of  the  summ 
frequency  (cosT=2n77)  and  'mf  is  the  number  of  cycles  of  the 
difference  frequency  (cocjT  =  2m77)  .  The  normalized  frequency 
errors  Aw  i  and  Aw  q  are  modeled  as  statisticlly 
independent  random  variables  uniformly  distributed  over 
[-0.1,0.1],  The  results  of  required  SNR  for  a  given  Pe  at 
different  values  of  JSR  are  summarized  in  Table  2. 


TABLE  2 
SNR  REQUIREMENTS  (P£  =  10" 6  FSK) 

JSR 

DETf  JAMM 

JAMM  WITH 
TIME  ERR 

JAMM  WITH 
FREQ'  ERR 

JSR=0.0 
JSR=0.1 

13.0  dB 
18.0  dB 

13.0  dB 
17.0  dB 

13.0  dB 
17.8  dB 

C.   INCOHERENT  RECEIVERS 

In  Chapter  3,  section  C,  it  was  shown  that  the  jammer 
model  utilized  which  includes  timing  errors  affects  the 
performance  of  incoherent  receivers  in  exactly  the  same  way 
as  a  jammer  with  no  timing  errors.  That  is,  the  timing  error 
associated  with  the  jammer  does  not  affect  the  resulting 
performance   of   the   incoherent   receiver.   Therefore   the 
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performance  differences  associated  with  the  jamming 
waveforms  with  frequency  errors  and  without  will  be  analyzed 
here . 

1.  ASK 

The  performance  of  the  incoherent  ASK  receiver  is 
graphically  displayed  in  Figure  (5.12).  These  plots 
correspond  to  numerical  evaluation  of  Equation  (3.9)  under 
the  assumption  of  equally  likely  hypotheses.  This  figure 
clearly  shows  the  'break  point'  effect  at  JSR  of  0.25.  It 
can  be  observed  that  16  dB  of  SNR  is  required  in  order  to 
obtain  a  Pe  of  10  when  no  jammer  is  present  (that  is 
JSR=0),  while  it  takes  23.5  dB  of  SNR  to  obtain  the  same  Pe 
for  a  JSR  value  of  0.1. 

Figure  (5.13)  displays  the  performance  of  the  same 
ASK  receiver  except  that  now  the  jammer  exhibits  frequency 
errors.  Figure  (5.13)  was  obtained  under  the  assumption  that 
the  number  of  waveform  cycles  per  observation  period  is  5 
and  the  normalized  frequency  error  A^m  can  again  be  modeled 
as  uniformly  distributed  random  variable  over  [-0.1,0.1]. 

Equation  (3.35)  was  used  to  generate  the  plots  of 
Figure  (5.13).  It  can  be  observed  that  at  JSR=0.1  in  order 
to  obtain  Pe  of  10"5'  the  value  of  SNR  required  is  23.0  dB . 
The  key  results  for  this  case  are  summarized  in  Table  3. 

2.  FSK 

Figure  (5.14)  and  Figure  (5.15)  correspond  to 
graphical  results  on  the  performance  of  the  coherent  FSK 
receiver.  When  the  jamming  waveform  is  modeled  as 
deterministic,  Pe  has  been  evaluated  using  Equation  (3.24) 
and  graphically  displayed  in  Figure  (5.14).  Here,  a  'break 
point'  occurs  at  a  JSR  value  somewhere  between  0.5  and  1.0 
as  shown  on  the  figure.  It  can  be  noted  that  13.5  dB  of  SNR 
is  required  in  order  to  obtain  a  Pe  of  10  for  a  JSR  value 
of  0.0,  while  the  same  Pe  is  obtained  by  increasing  the  SNR 
to  16.5  dB  for  a  JSR  value  of  0.1. 
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TABLE  3 


-5 


SNR  REQUIREMENTS  (PE  =  10 ~D   ASK) 


JSR 

DET'  JAMM 

JAMM  WITH 

JAMM  WITH 

TIME  ERR 

FREQf  ERR 

JSR=  0.0 

16.0  dB 

16.0  dB 

16.0  dB 

JSR=  0.1 

23.5  dB 

23.5  dB 

23.0  dB 

JSR=0.25 

unable 

unable 

large 

Figure  (5.15)  displays  Pe  plots  corresponding  to  the 
same  receiver  under  the  assumption  that  tha  jammer  contains 
frequency  errors.  In  this  case,  16.8  dB  of  SNR  is  needed  in 
order  to  obtain  a  Pe  of  10"5  at  a  JSR  value  of  0.1.  The  key 
results  are  summarized  in  Table  4. 


TABLE  4 
SNR  REQUIREMENTS  (P£  =  10 "5  ASK) 


JSR 

DET?  JAMM 

JAMM  WITH 

JAMM  WITH 

TIME  ERR 

FREQ'  ERR 

JSR=0.0 

13.5  dB 

13.5  dB 

13.5  dB 

JSR=0.1 

16.5  dB 

16.5  dB 

16.8  dB 

JSR=0.3 

23.5  dB   123.5  dB 

1 .  i 

22.5  dB 

These  results  were  obtained  under  the  assumption  that  the 
normalized  frequency  errors  Aa>jyj  i  and  Aco^j  q  are 
statistically  independent  random  variables  uniformly 
distributed  over  [-0.1,0.1]  and  n=6,  m=4  while  Tn'  is  the 
number  of  cycles  of  summ  frequency  per  observed  period  and 
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'm'  is  the  number  of  cycles  of  the  difference  frequency  per 
observed  period. 
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V.  CONCLUSIONS 

In  this  thesis,  the  effect  of  timing  errors  or  frequency- 
errors  associated  with  jamming  waveforms  used  to  degrade  the 
performance  of  coherent  as  well  as  incoherent  receivers  has 
been  analyzed.  It  has  been  assumed  that  these  receivers  are 
designed  to  operate  on  a  signal  plus  white  Gaussian  noise 
only  environment . 

The  analysis  of  the  effectiveness  of  jammers  that 
exhibit  timing  and  frequency  errors  was  undertaken  by 
evaluating  the  receiver  probability  of  error  in  the  presence 
of  jamming,  and  comparing  the  results  to  those  obtained  when 
the  jamming  had  no  timing  or  frequency  errors.  Thus, 
receiver  error  probabilities  as  a  function  jammer- to- signal 
ratio  (JSR)  and  signal-to-noise  ratio  (SNR)  were  evaluated 
and  determinations  were  made  on  how  much  larger  JSR  values 
were  required  as  a  result  of  timing  or  frequency  errors,  in 
order  to  produce  a  given  Pe  at  some  fixed  SNR  value. 

For  the  coherent  receivers  analyzed,  the  effect  of  the 
timing  error  associated  with  the  jammer  was  studied  and  it 
was  concluded  that  the  jammer's  lack  of  sychronism  with  the 
digital  data  causes  a  small  decrease  in  jammer 
effectiveness.  Typically  1  or  2  dB  more  JSR  was  needed  in 
order  to  overcome  timing  errors.  Due  to  the  lack  of  correct 
knowledge  of  the  frequencies  of  communications,  the  jammer 
required  typically  1  or  2  dB  more  power  in  order  to 
compensate  for  frequency  errors. 

For  the  incoherent  receivers  analyzed,  it  was  determined 
that  the  timing  errors  associated  with  the  jamming  waveform 
have  no  effect  on  the  performance  of  the  receiver.  However, 
the  jammer's  incorrect  knowledge  of  the  communication 
system's  operating  frequencies  results  in  a  loss  of  jammer 
effectiveness   to   the   point   that   jammer   power   must   be 
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incresed  by  about  1  to  2  dB  in  order  to  overcome  the  lack 
of  precise  knowledge  of  the  communication  system's  operating 
frequencies . 
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Figure  5.4    Incoherent  Matched  Filter  Receiver 
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Performance  of  the  PSK  Coherent  Receiver 
with  Deterministic  Jamming 
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Performance  of  the  PSK  Coherent  Receiver 
with  Timing  Error  Jamming 
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Figure  5.8    Performance  of  the  PSK  Coherent  Receiver 
with  Frequency  Error  Jamming 
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Performance  of  the  FSK  Coherent  Receiver 
with  Deterministic  Jamming 
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Performance  of  the  FSK  Coherent  Receiver 
with  Timing  Error  Jamming 
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Figure    5.11         Performance    of    the    FSK    Coherent    RECEIVER 
with   Frequency    Error   Jamming 
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Figure    5. 12 


Performance    of    the   ASK    Incoherent    Receiver 
with   Deterministic    Jamming 
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Figure  5.13 


Performance  of  the  ASK  Incoherent  Receiver 
with  Frequency  Error  Jamming 
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Figure  5 . 14 


Performance  of  the  FSK  Incoherent  Receiver 
with  Deterministic  Jamming 


84 


FSK  (FREQ'  ERROR) 


ioo  20.0 

SNR(DB) 

INCOHERENT  RECOVER 


4O0 


Figure  5.15    Performance  of  the  FSK  Incoherent  Receiver 
with  Frequency  Error  Jamming 
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